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Abstract 

The results of the recent experiments on the reaction t:~p — > vr'^vr^n performed 
at KEK, BNL, IHEP, and CERN are analyzed in detail. For the / = vrvr 5 wave 
phase shift 5q and inelasticity t/q a new set of data is obtained. Difficulties emerging 
when using the physical solutions for the S and D wave amplitudes extracted 
with the partial wave analyses are discussed. Attention is drawn to the fact that, for 
the TT^TT^ invariant mass, m, above 1 GeV, the other solutions, in principle, are found 
to be more preferred. For clarifying the situation and further studying the /o(980) 
resonance thorough experimental investigations of the reaction ■K~p -K^ir^n in the 
m region near the KK threshold are required. 
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I. INTRODUCTION 



The reactions ttN — > 'K'kN so far are a major source of information on the processes 
TTTT — > TTTT. At high energies and small values of the momentum-transfer-squared from the 
incident vr to the outgoing vrvr system, < — t < 0.2 GeV^, the reactions vrA^ nnN are 
dominated by the one-pion exchange (OPE) mechanism. In treating the data on these 
reactions the partial wave analysis method is used. As a rule, a few possible solutions 
for the partial wave amplitudes of the final tttt system are obtained. In some cases, 
the preferred solution is selected from the additional physical arguments. Generally, 
to obtain the reliable and unambiguous results in a wide region of m, high statistics, 
polarized targets, and precise measurements of the ttN rniN reaction cross section 
at different energies are needed. The detailed reviews and comprehensive discussions of 
the experimental results on the reactions ttN — > ttttN and tttt scattering in the region 
2m^ < m < 2 GeV available by the early 1999 have been presented in Refs. [1,2]. 

In this work we analyze the recent data on the intensities and relative phase of the 
S and D partial waves of the tt'^tt'^ system produced in the reaction n'p n^TT^n. The 
data have been obtained in the experiments with incident 7i~ energies of 8.9, 18.3, 38, 
and 100 GeV performed at KEK [3], BNL [4], IHEP [5], and CERN [6], respectively. 
Our main goal is to obtain information on the tttt S wave phase shift 5q and inelasticity 
r/o in the channel with isospin / = that would be complementary to the previous 
"canonical" data extracted from the 17.2 GeV experiments on the reactions n^p 'n:^'K~n 
[7-12]. We especially emphasize a strong likeness of the physical solutions selected in all 
four experiments on tt'^tt^ production and also the common difficulties that emerge when 
interpreting these solutions and in their comparison with the 7r"'"7r~ data. It turns out, 
in particular, that some of the physical solutions found lead to considerable violations 
of the unitarity condition for the tttt scattering amplitude in question. In addition, we 
conclude that the data on 7r°7r° production are indicative of a noticeably smaller value of 
the /2(1270) TTTT decay branching ratio in comparison with the Particle Data Group 
(PDG) data [13]. In connection with a considerable interest in the fight scalar meson sector 
(see for reviews Refs. [1,2,13-15]), we suggest to perform especially careful measurements 
of the reaction tt~p ir^TT^n in the m region from 0.9 to 1.1 GeV, i.e., near KK threshold. 
This would allow the /o(980) coupling constant to the KK channel to be determined more 
reliably and also to resolve the long-standing question [16] of a possible ambiguity in the 
behavior of the phase shift 5° above the KK threshold. 

The paper is organized as follows. In Sec. II, the KEK results [3] are analyzed. In Rcf. 
[3] the data on the phase shift 5q have been obtained in the m region from 0.36 to 1 GeV. 
The 5q values found by us in other way in the interval 0.68 < m < 1 GeV agree with the 
KEK data [3] within experimental uncertainties. We also present new results for 5^ and 77° 
in the region 1 <m < 1.64 GeV. In Sec. Ill, the extrapolation of the S and D wave mass 
distributions obtained in the BNL experiment [4] from the physical region of the reaction 
7r~p n^n^n to the pion pole (t = m^) is performed. Considering the different solutions 
found in Ref. [4] for these distributions we obtain a few sets of the values of 6q and t^q 
in the m region from 0.32 to 1.52 GeV. The GAMS results on the reaction n^p — > 7r°7r°n 
[5,6] are discussed in Sec. IV. Here we also summarize all difficulties encountered while 
analyzing the tt^tt'' data [3-6]. In Sec. V, we formulate briefly a few concrete suggestions 
for further studying the reaction TT~p n^n^n which, as one can hope, will be used to 
clarify the experimental situation. 
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II. ANALYSIS OF THE KEK DATA 



In the KEK experiment [3] , the data on the intensities and the relative phase of the S 
and D partial waves for the reaction vr+vr" tt^tt^ have been obtained in the m interval 
from 0.36 to 1.64 GeV. They have been extracted from the n~p — > 7r°7r°n data by using 
the linear Chew-Low extrapolation and partial wave analysis. Because the absolute 7r°7r° 
production cross section has not been determined in the experiment, the S and D wave 
intensities, |^sp and |^dP, have initially been presented in arbitrary (identical) units 
[3]. Any alternative solution for jAgp and jA^ip has not been discussed in Ref. [3]. The 
S and D wave intensities arc related to the phase shifts 6q and 62 and inelasticities tJq 
and in a conventional way: \As\'^ ~ {a^ — CqP, where Oq = (?7q exp(2i5Q) — l)/2i, and 
I^Dp ~ \a°2 - a^P, where = (77I exp(2i5|) - l)/2i. To find (5° below the KK threshold, 
it was assumed [3] that in this region ?7q = ?7q = 1, and consequently |^sP ~ sin^((5Q — 6q). 
As is well known from a large number of previous experiments, the phase shift 6q smoothly 
goes through 90° in the region 0.7 < m < 0.9 GeV and the phase shift 6q is negative, 
smooth, and small (see, for example, [2,7,8,17]). Therefore, to extract the phase shift 
difference 5q — 5^ from the unnormalized data, the following normalization condition was 
accepted in Ref. [3]: the maximum value of l^sp is equal to 1. The KEK data for the S 
and D partial wave intensities normalized in this way ara shown in Fig. 1, together with 
the data on the relative phase S = 4>s — 4>d between the amplitudes As = \As \ exp {i4>s) 
and Af) = \Aj:, \ exp {i(j)D)- The values of the I — 2 nn S wave phase shift 6^ used in Ref. 
[3] were given by the parametrization 6q — —0.87q [with q — (m^/4 — m^)^/^ in GeV and 
Sq in radians], and in such a way the data for Sq were obtained in the m region from 0.36 
to 1 GeV. In the following, for 6q we shall use the fit to the data from Refs. [17,18] which 
is shown in Fig. 2. ^ Using the fit and the data for [Asp shown in Fig. la we have also 
determined the values of Sq for m < 1 GeV. They are plotted in Fig. Id. The resulting 
values are in excellent agreement with those obtained in Ref. [3]. 

We now determine 5" and simultaneously by using the available data on the relative 
phase S and the intensity {Agl"^ (see Figs. Ic and Id) in the m region from 0.68 to 1.64 
GeV. In order to estimate the phase 4>d 'Vfe neglect the tiny amplitude [17,18] (which 
is quite reasonable because the experimental errors of as is seen from Fig. lb, 

are not too small) and assume that the D wave amplitude is dominated by the /2(1270) 
resonance contribution and can be written in the form 

— — 2 2 '■ F' ^' 

mj^ — m'^ — imf^V 

were F = {mjJm)T f,{q/ qf,fD{qM / D{qRf,), D{x) = 9 + 3^^ + x^ qf, = {m)jA- 
m^)-*^/^, Rf^ is the interaction radius, and rrif^, Fj^, and -B/jtttt are the mass, width, and 
TTTT decay branching ratio of the /2(1270). The fitted curve on Fig. lb corresponds to the 
following values of the /2(1270) resonance parameters: 

mj, = 1.283 ±0.008 GeV, F/^ = 0.170 ±0.014 GeV, 

7?/, = 3.59 ± 0.71 GeV-\ S/^^^ = 0.760 ± 0.034 . (2) 

-"^This fit was obtained with the parametrization 6q = —aq/{l + hm? + cm^ + dm^), where a, b, c , 
and d are fitted parameters. In addition, in lacking the reliable data on the deviation of tjq from 1, we 
set ?7g = 1 for all considered values of m. When more or less detailed data on ?7q arc available, it will be 
interesting to take into account possible inelastic effects. It is not unreasonable to think that such effects 
may appear in the I = 2 rnr channel only above the nominal pp threshold (1.54 GeV), but not above the 
KK threshold as in the I = ww channel. 
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Thus, the phase is defined by that of the Breit-Wigner amphtude (1). To express the 
parameters 6q and rj^ in terms of the known values S, Sq , and it is convenient 

to represent the amphtude As in the form (see footnote 1) 

As = e^^^o ^ e'<As = e^(^^o+<^)|I^| , (3) 

where (f) is the phase of the amphtude As- The distinctive feature of the amphtude As 
is that in its Argand diagram the relations between 5q — 5q, (p, 7)^, and \As\ formally 
look like the relations between the corresponding parameters of any unitary partial wave 
amplitude with definite isospin /; for example, the phase is confined within the range 
from 0° to 180° because lm{As) > 0. Thus, we have 

<P^5-25l + <f>D, 77° = Vl- 41^51 sin + 4|A5|2, (4) 



• ofxo x-2\ 21^51 COS0 2^ l-2\As\sm(f) 

sm2(5^-do) = , cos2(5^-5o) = ^ . (5) 

Vo Vo 

Since the interference between the S and D partial waves is defined by the product 
1^451 174^| cos 5 and cos 5 determines S only up to the sign, two solutions always exist: the 
solution with 5 > and the other one with S < 0. Moreover, if cos S is close to 1 (and 
|5| 0) in some region of m then in this region a transition from one solution to the other 
one is possible. The KEK data [3] presented in Fig. Ic show that the phase 5 changes 
most rapidly near the KK threshold (which is one of the evident manifestation of the 
/o(980) resonance) and that just near 1 GeV cos 5 ^ 1. Thus, in principle, we have four 
possible variants: (i) 5 > for m < 1 GeV and 5 < for m > 1 GcV, (ii) S > for 
all m, (iii) 6 < for all m, and (iv) 6 < for m < 1 GeV and 6 > for m > 1 GeV. 
However, variants (iii) and (iv) with 5 < for m < 1 GeV can be rejected at once. Really, 
estimating 5 for m < 1 GeV by using the relation S — Sq + 6q — 4>d, one can easy verify 
that 5 must be positive in this region with the conventional definition of the signs of the 
phase shifts (5°, 5q, and [see Figs. Id, 2, and Eq. (1)]. So, we shall consider only 
variants (i) and (ii). 

Figure 3 shows the values of 6q and extracted from the KEK data (sec Figs, la, 
lb, and Ic) in the region 0.68 < m < 1.64 GeV by using the Eqs. (4) and (5) for the 
two above mentioned variants of the 5 phase behavior. The values of 5q in the region 
0.36 < m < 1 GeV obtained above from the data on \As\'^ with r]^ = 1 (sec Fig. Id) are 
also shown in Figs. 3a and 3b for comparison and completeness. As is seen, for example, 
from Fig. 3a, the sets of the 6^ values found in the region 0.68 < m < 1 GeV by two 
different ways are in quite reasonable agreement with each other. With obtaining 6q and 
rjQ in the general case, the Argand diagram of the amplitude ^4^ was built for each variant. 
After this the values of 2(5° — Sq) obtained from Eqs. (4) and (5) were finally defined 
by the requirement that those of 6q be smoothly connected as a function of m. That the 
strong violation of unitarity takes place in variant (ii) for m > 1.16 GeV (see Fig. 3d) 
can be easily understood from the relation (p = 6 — 26q + 0d [see Eq. (4)]. The fact is 
that the values of in this case fall into the range from 180° to 360°, which is forbidden 
as (f) is the phase of the formally unitary amplitude As- Furthermore, in variant (ii), the 
phase sift Sq for m > 1 GeV (see Fig. 3c) is in rather poor agreement with the tt+tt" 
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production data [7-12] according to which, for example, at m ^ 1.3 GeV Sq has to be 
close to 270°. Thus, variant (ii) with 5 > for all m can be rejected. As for variant (i), 
there are a set of specific features which, to our knowledge, are missing from the n'^n' 
data [7-12]. As is seen from Figs. 3a and 3b, in this case we have noticeable differences 
of tJq from unity for m < 1 GcV, its approximate equality to 1 for 1 < m < 1.12 GeV, 
violation of unitarity near 1.2 GcV, and sharp jumps of 6q and rjQ with further increasing 
m. There is little doubt that these features are artefacts of the partial wave analysis of 
the 7r~p — > 7r°7r°n data. 

Another difficulty is that the accepted normalization for \As\'^ leads to -B/jtttt = 0.760± 
0.034 [see Eq. (2)], while according the PDG data [13] Bf,^^^ = 0.847±g:g?|. These two 
values cling to one another only by their double errors. Hence, in principle, one may 
conclude that the tt^tt^ data [3] indicate that the absolute cross section of the /2(1270) 
resonance formation through the OPE mechanism in the reaction 7r~p — > 7r°7r°n can turn 
out to be approximately 20% smaller, at least at m ~ m/j, than that expected from the 
PDG data [13]. Alternatively, the KEK data [3] might be normalized with use of the 
known value max|A/)p = (1 +7^0)^/4 = Bj^^^ with Bj.^^^^^ from Ref. [13]. However, in 
this case, the resulting values of \As\'^ in the most interesting region of the lightest scalar 
resonance cr(600) [3,13,19] would be approximately 25% higher than the unitarity limit 
for 1'- 

We shall see in the next Sections that the other experimental data on the reaction 
Tr~p — > 7r°7r°n lead to very similar difficulties. 

III. ANALYSIS OF THE BNL DATA 

In the BNL experiment [4], the high statistics on the reaction 7r~p Ti^ir'^n (about 
8.5 X 10^ events) has been accumulated and the detailed partial wave analysis of the 7r°7r° 
angular distributions has been performed. This analysis has been done for ten sequential 
intervals in —t covering the region < —t < 1.5 GeV^ and over the m range from 0.32 
to 2.2 GeV scanned with the 0.04 GeV -wide step. As a result, two solutions for the 
unnormalized intensities of the S and Dq partial waves and four ones (because of a sign 
ambiguity) for their relative phase have been obtained. The above quantities were denoted 
in Ref. [4] by j^j^, j-Doj^, and ips-Do, respectively; in so doing, Dq denotes the D wave 
with Lz — 0, where is a projection of the 7r°7r° relative orbital angular momentum 
on the z-axis in the Gottfried- Jackson reference frame [4]. In the following, we shall use 
these notations, too. One of the solutions for the S and Dq wave intensities, which is 
characterized by a large magnitude of 15*1^ and a small one of |-Do|^ for ni < 1 GcV, and 
which is smoothly continued to the higher mass region, where the Dq wave is dominated 
by the /2(1270) resonance contribution, has been selected in Ref. [4] as the physical 
solution. Together with the intensities \S\^ and \Dq\'^, the physical solution also includes 
two corresponding sets of the ^s-Dq phase values which differ only in sign. Note that 
the other solution intersects with the physical one aX m \ GeV. We agree with the 
physical arguments given in Ref. [4] based on which the other solution can be rejected in 
the region m < 1 GeV. However, for m > 1 GeV we shell analyze the two solutions and 
also the cases with transitions of the phase ^ps-Do from one solution to the other one. 

For the analysis we take the BNL data [4] on 15*1^, |-Do|^, and <fs-Do pertaining to five 
intervals of -t, 0.01 < -t < 0.03 GeV^, 0.03 < -t < 0.06 GeV^, 0.06 < -t < 0.1 GeV^, 
0.1 < -t < 0.15 GeV^, and 0.15 < -t < 0.2 GeV^, and to the region 0.32 < m < 1.6 
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GeV. Note that the data on (ps-Oo available only for m > 0.8 GeV. To obtain the 
values of the quantities l^ijp, and 5 (see Sec. II) as functions of m characterizing 

the reaction tt'^tt" — > 7r°7r° on the mass shell, we parametrize the t dependence of \S\'^, 
|DoP, and (fis-Do by means the following expressions 

ic^|2 |2 -^exp[&g(t-m^)] 2 .^\a |2 -^exp[6z)o(t-m^)] 

P = \^S\ 77 , \^o\ =0 — \Ad\ 7- ^ , [b) 

q [t- m^y q [t - m^y 

<Ps-Do + a {t/ml - 1) (7) 

and, in each 0.04 GeV mass bin, thus determine the unnormalized intensities 1^5^ and 
the phase 6 and also the slopes bs, &Do, and a by fitting to the BNL data on the t 
and m distributions by the formulae (6) and (7). ^ In so doing, for jS'p, |-DoP, and <ps-Do 
in each —t bin we take into account the physical solution for m < 1 GeV and the physical 
and other ones for m > 1 GeV. Unfortunately, the absolute value of the 7r~p — > 7r°7r°n 
reaction cross section has not been determined in the BNL experiment [4]. Therefore, to 
normalize the extrapolated intensities jA^p and jA/jp we proceed in the same way as in 
Sec. II. The extrapolated and normalized data corresponding to the physical and other 
solutions are plotted in Figs. 4a, 4c, and 4e with solid and open symbols, respectively. 
It is interesting to note that as a result of the extrapolation two branches of the (fs-Do 
phase values for the other solution (i.e., the branch with (fis-Do > for all m and that 
with <fs-Do < for all m) interweave with each other in the region m > 1.24 GeV (sec 
Fig. 4e) in such a way that there arise two new branches of the extrapolated phase 6, 
which are characterized by a smooth dependence on m and which, for example, can be 
considered either intersecting or osculating near 1.26 GeV. 

As in Sec. II, we begin with the determination of the phase shift 6^ for m < 1 GeV 
from the data on (see Fig. 4a) assuming that 77° = 1 in this region. The resulting 

phase shift values are shown in Fig. 5 by open circles. Note that two points in the region 
m ~ rriK disturbed by the Kg — >• 7r°7r° events [4] arc omitted. Then, heaving the data 
on I^Dp, and 5 for m > 0.8 GeV (see Fig. 4), we determine the values of Sq 

and with use of the general formulae (3) and (4), and also Eq. (1). The results for 
the previously selected solutions among all the possible ones, that are shown in Fig. 4, 
are plotted in Fig. 5 with solid circles. Strictly speaking, the selection is reduced to 
rejection of the physical solution with 5 < for m < 1 GeV (see Fig. 4e), since a simple 



^Such two-parametric fits to the off-shall partial wave intensities were widely used in the literature to 
obtain the suitably extrapolated data (see, for example, Refs. [9,17,20,21]. However, the determination of 
the phase S with use of the direct extrapolation of the data on the phase ips-Da [sec Eq. (7)] may provoke 
a question. If the data on the S — Dq interference contribution, as such, had been presented in Ref. [4], 
the problem would not have arisen. The fit to such data to the function —2taexp[b{t — m^)]/(f — m^)^ 
analogous to those in Eq. (6) and the identification of the fitted parameter a with \/5('7?-^/'i')|^s| cos S 
would allow I (5 1 to be determined in the proper way. Because such data are not available, the indirect test 
of the results obtained with Eq. (7) was carried out. Using the data [4] on jSI^, |£>o|^) and fs-Do we 
constructed the quantity 215*1 |_Do| cosips-Do and foimd with the above extrapolation the on-shcU S — D 
interference contribution. Then, knowing independently \As\'^ and we determined (5 as a function 

of m. The 5 phase values obtained in the two ways are in very close agreement with each other. Certainly, 
owing to the forced double rccoimting of the errors of the input data with the indirect test, the errors of 
S turn out to be larger then those obtained from the fit by Eq. (7). On the other hand, when the values 
of 6 are determined by using Eq. (7) their errors practically are not differ from the errors of the input 
data for tps-Do- AH the aforesaid allowed us to prefer the determination of the phase 6 with use of Eq. 
(7). 
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estimate 6q = 6 — 6q + (pD (see Sec. II) yields Sq ~ —(25 ^ 40)° for this solution in the 
region of 0.8 — 1 GeV, which is, certainly, unsatisfactory. In its turn. Figs. 5a and 5b 
show that the physical solution with 5 > for all m can be also rejected due to strong 
violation of the unitarity condition for m > 1.2 GeV. Figures 5c and 5d correspond to 
the physical solution for l^spj l^-Dpj and S with the transition of the phase 5 at m ~ 1 
GeV from the branch with 5 > to that with d < (sec Fig. 4c). Such a physical 
solution consists with unitarity but corresponds to the weak coupling between the tttt and 
KK channels near the KK threshold. Indeed, for this solution rj^ is close to unity in the 
region 1 < m < 1.15 GeV. However, the latter disagrees with the data obtained from the 
reactions n^p 7r~^iT~n, 7r~^p 7r"'"7r~ A^"*", and nN KK{N, A) (see, for example, 
Refs. [7,9,10,16,22,23]). Figures 5c and 5f correspond to the combination of the physical 
solution with 5 > for m < 1 GeV and the other one for m > 1 GeV with 6 > in 
the region 1 < m < 1.28 and 5 < in the region 1.28 < m < 1.52 GeV (see Figs. 4a 
and 4e). Finally, Figs. 5g and 5h correspond to the similar combination of the physical 
solution and the other one for which S < for m from 1 to 1.52 GeV (see also Fig. 4e). 
Certainly, there are two more variants which differ from the last two ones only by the 
sign of 5 for m > 1.28 GeV (see Fig. 4e). These variants lead, however, to appreciable 
violations of the unitarity condition for m > 1.32 GeV, and therefore, are of little interest. 
Thus, one can conclude that just the variant presented in Figs. 5e and 5f is, in many 
respects, in qualitative agreement with the results of the previous partial wave analyses 
of the TT+TT^ data [1,7,9,11]. As indicated above, this variant corresponds to the positive 
relative phase 6 = 4>s~4>d up to the /2(1270) resonance and the negative one above it. An 
important point is that such a behavior of 5 as a function of m is strongly confirmed by 
the pioneering data from the polarized target experiment on the reaction Tr~p — > 7r'*"7r~n 
at 17.2 GeV [11]. 

Comparing Fig. 5 with Fig. 3 we just note that the BNL data lead to the obviously 
higher values of the phase shift 5^ for m < 0.5 GcV than the KEK data. 

In extracting the information on Sq and tIq, the phase (pD has been defined by fitting 
to the data on \Ad\'^ (see Fig. 4c) with use of Eq. (1). The parameters of the /2(1270) 
were found to be (see also the curves on Fig. 4c): for the physical solution, 

= 1.279 ± 0.002 GeV , F/^ = 0.205 ± 0.005 GeV , 

Rf, = 3.96 ± 0.24 GeV-^ , Bf^^^ = 0.697 ± 0.008 (8) 

and, for the other solution, 

= 1.281 ±0.002 GeV, F/^ = 0.211 ±0.005 GeV, 

^ 4.65 ± 0.33 GeV~^ , Bf^^^ = 0.712 ± 0.007 . (9) 

Thus, the BNL data indicate that the branching ratio Bf^j^j^ can amount to approximately 
84% of the PDG value [13]. Possible consequences of a similar discrepancy has been 
already discussed in connection with the KEK data at the end of Sec. II (recall that the 
relevant ratio for the KEK data has been found to be approximately 90%). 

IV. DISCUSSION OF THE GAMS DATA 

The highest statistics on the reaction n~p — > n^n^n was accumulated by the GAMS 
Collaboration in two experiments at 38 GeV [5] and 100 GeV [6]. However, the small —t 
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region from to 0.2 GeV^ has been examined in the works [5,6] very sparingly. But the 
data averaged over the small —t region have been presented for l^'p and (fis-Do iii R-^f. [5] 
and for l^'p, |-DoP, and (ps-Do in Ref. [6]. Such "global" data, of course, do not permit 
to perform a proper extrapolation of the mass distributions measured from the physical 
region to the pion pole. Nevertheless, we discuss some typical features of the GAMS data. 
The physical solution for \S\'^ and <fs-Do and the other solution only for [S'p have been 
presented in Ref. [5] in the region 0.8 < m < 1.6 GeV. The (fs-Do phase was found to be 
positive in the full mass range [5] (about existing the ambiguous solution with (ps-Do < 
the readers, probably, have to guess by themselves). In general, the available GAMS data 
[5] are very similar to the corresponding BNL data [4]. For example, in the case of the 
physical solution, \S\'^ and ^s-Dq from Ref. [5] behave as functions of m in the same 
way as the extrapolated quantities \As\'^ and 5 shown in Figs. 4a and 4e by solid circles. 
However, it is such a physical solution for \As\'^ and 5 (with 5 > for all m) that leads to 
strong violation of the unitarity condition for m > 1.2 GeV (see Fig. 5b). In analyzing 
the 7r°7r° system produced in the reaction 7r~p — > 7r°7r°n at 100 GeV, the only solution 
for IS'p, l-DoP) and Vs-Dq has been selected and presented by the GAMS Collaboration 
in Ref. [6]. Unfortunately, this unique solution is very close to the above physical one 
obtained in the GAMS 38 GeV 7r~p — > 7r°7r°n experiment [5]. 

It is of first importance that the GAMS GoUaboration measured the absolute cross 
section of the /2(1270) resonance formation in the Dq wave in the region < — t < 0.2 
GeV2. According to Ref. [24], at 38 GcV, aooiTr-p /2(1270)n ~> = 2.3 ± 0.2 

/A). More recently, this value was used, in particular, to normalized the 100 GeV data [6]. 
Although the cross section value obtained is approximately 1.5-2 times greater than in a 
set of previous n^n^ production experiments [24,25], nevertheless, it is 1.57 times smaller 
than the estimate based on the OPE model (it is an old story with the experimental 
underestimation of the 7r~p /2(1270)n — > TT^n^n reaction cross section all details of 
which can be found in Ref. [25]). By using this model with the PDG values of m/j, Tf^, 
and -B/jTTTr [13] we estimate 

(^Do(t^~P /2(1270)n ^ 7r°7r°n) ^ a^^^in'p /2(1270)n ^ 7r\%) ^ 

-o.2GeV^ 

where P^- = 38 GeV, ^ ^ 2x14.3, andft/^ ^ 7.5 GeV-2+2x0.8 GeV'^ ln(38/18.3) 
~ 8.68 GeV^^ (in estimating the slope bj^, its Regge energy dependence and the results 
for the slope b^^^ in the /2(1270) mass region presented in Fig. 4d have been taken into 
account). Note that this estimate is in good agreement with the result of the extrapolation 
of the available data on the reaction n^p /2(1270)n — > 7r+7r~n at 17.2 GeV [8], 100 
GeV and 175 GeV [26] to the GAMS energy (see Ref. [25] for details). Thus, the GAMS 
data [24] indicate that the value of P/jtttt can amount to about 80% of that given by the 
PDG [13]. 

We now summarize briefly the common difficulties encountered in analyzing the data 
obtained in four recent experiments on the reaction 7r~p — > tt^tt'^ti. First, the physical 
solutions selected by using the partial wave analyses of the tt^tt^ production data lead to 
the values of Sq and tJq which are incompatible with the known results obtained from the 
7r+7r~ data, at least for m > 1 GeV. Some of these solution lead to strong violations of 
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the unitarity condition. On the other hand, among the other solutions one can point out, 
in principle, the more preferred ones. Secondly, it is astonishing that the data of the four 
recent experiments on 7r°7r° production include the indications that the value of -B/jTrTr 
can be distinctly smaller than the currently accepted one. This difficulty is rather serious 
and highly interesting. Let us remind that the tttt production experiments on unpolarized 
targets, in particular, those under discussion here, do not permit the contributions of the 
TT and ai exchange mechanisms to be separated in principle, even with huge statistics, 
because these contributions to the unpolarized cross section are incoherent [27]; in other 
words, there is no model-independent way to do this. Therefore, by our opinion, the 
difficulty with Bf^jm may present itself a further evidence that the partial wave analyses 
of the unpolarized data allow to determine the intensities and the relative phases of the S, 
D, ... Tin partial waves only approximately, with any extrapolation method. "The degree 
of proximity" is associated with the relative magnitude of the non-leading ai exchange 
contribution. With high statistical accuracy of the unpolarized data the presence of the 
tti exchange mechanism can manifest itself in the events responsible for \S\'^ just in the 
form of the above difficulty. In fact, this statement follows naturally from the analysis of 
the unnormalized KEK and BNL data (sec the end of Sec. II and also Sec. V for details). 
As for the GAMS data [23] , they appear to point merely to the general problem involving 
the accurate measurement of the Tr~p — > 7r°7r°n reaction cross section. More discussions 
both of the additional assumptions needed in analyzing the unpolarized data and of the 
oi exchange contribution can be found in the works [1,2,11,12,27-29]. 

V. CONCLUSION 

Using the most simple way we have extracted the values of the I = tttt S wave 
phase shift 5q and inelasticity tJq from the current data on the reaction 7r~p 7r°7r°n. 
It seems clear that a new set of precise experiments on this reaction is needed both for 
the more precise definition of the tt^tt" production mechanism and for obtaining the more 
detailed information on tttt scattering and light scalar resonances in the tttt channel. Let 
us formulate in this connection several concrete suggestions, leaving aside the general wish 
to investigate the reaction 7r~p 7r°7r°n on the polarized target. 

1) Detailed data on the m and t distributions for the tt^tt^ S and Dq partial waves, 
especially in the region < — t < 0.2 GeV^ where the OPE mechanism dominates, and 
measurements of the absolute value of the 7r"p 7r°7r*^n reaction cross section at different 
energies, for example, at KEK, BNL, IHEP, and CERN, would be highly desirable. The 
relative accuracy of new measurements must be comparable with (or better of) that given 
by the PDG [13] for Bj^^^. This would allow to perform the accurate description of the 
/2(1270) formation differential cross section within the OPE model and to test how well 
the S wave vr^vr" production cross section at its absolute maximum (which is located in 
the region 0.6 < m < 0.8 GeV) agrees with the OPE model prediction under the standard 
normahzation condition according to which 1^5^ = 1 (i.e., Sq—Sq — 90° and tJq — tjq — 1) 
at the absolute maximum point. An excess of the experimental values over the model 
expectations would be a good evidence, obtained from the unpolarized target data, for 
the presence of the ai exchange contribution to the S wave tt^tt^ production cross section 
in the region of its absolute maximum. Alternatively, if the maximal experimental value 
of the S wave cross section turns out to be less than in the OPE model, then it will 
completely disturb of the existing ideas about the Sq phase shift for m < 1 GeV, which 
seems to be highly unlikely. 
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2) We suggest to perform in the low —t region the especially careful measurements of 
^0^0 production in the S wave for m from 0.9 to 1.1 GeV, i.e., in the region of the well 
known interference minimum in \S\'^ located near the KK threshold. This would allow to 
obtain the important additional information on the /o(980) resonance coupling constant 
to the KK channel, gj^KK^ and to resolve the long-standing question [16] concerning a 
possible ambiguity in the behavior of the phase shift 6q above the KK threshold which 
arises at gj^j^+j^-/4:7i > 47rm|- ^3.1 GeV^. Furthermore, the magnitude of the S wave 
intensity in the immediate region of the minimum (if it lies below the KK threshold) can 
be used to obtain a very strong upper limit on the ai exchange contribution at small —t 
in this region of m. 

3) As a rule, the assumption of phase coherence between the Dq and D_ amplitudes is 
one of those using to select the physical solution, see, for example, Refs. [4,6,9,30]. Here 
D_ denotes the D wave with |L^| = 1, in the Gottfried- Jackson reference frame, which is 
produced via the unnatural parity exchanges in the t channel of the reaction tcN ttttN. 
In this connection we would like to call attention to a new curious circumstance. According 
to the GAMS measurements [6], the ratio jD-lVl^oP in the /2(1270) mass region at 100 
GeV is half as large as that at 38 GeV. This fact may testify to compensation of the 
TiP and a2-P Regge cut contributions (P denotes the Pomeron exchange) to the D_ wave 
production amplitude with energy increasing, i.e., to violation of phase coherence. 

4) We also suggest to use in future the more suitable notation for the S — Dq in- 
terference contribution to be extracted in the unpolarized target experiments, instead of 
the commonly used simplified one of the form \S\\Do\ cos ips-Do- It includes the men- 
tion of the coherence factor (see, for example, Ref. [31]) and is more adequate to the 
measured quantity. Experimentally, the 5* and Dq wave intensities, 15*^ and jP'oP, and 
the S — Dq interference contribution, ,^ 15*1 |Do| cos are measured simultaneously. In 
fact, \S\ = [l^^p + |^aJT^^ l^ol = [I^OttP + |P'oaJT^^ ^ud the coherence factor ^ 
(0 < ^ < 1) and the phase if have the form 

= arctan[( E \Si\\Doi\sinifi)/ { \Si\\Doi\ cos (pi)] , 

where {Dq^^) and (DoaJ are the S (Dq) wave production amplitudes caused by 
the TT and ai exchange mechanisms, respectively (these amplitudes correspond to two 
independent configurations of the nucleon helicities in the reaction ttN — >■ irirN), and ipi 
is the relative phase between the amplitudes Si and D^i. Let us consider the case when 
the amplitude Doai is negligible. Then, denoting = ip,, hj ips-Do^ one can see that 
the real interference contribution differs from that presented with the simplified notation, 
|S'||Do| cos(y95_£)Q, by the coherence factor ^ = 1 / yT+~|S^^JY]5Vp . If we put ,^ = 1 for 
all m, we shall always deal with the effectively underestimated values of | cos(P)s-Dq\- 
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Figure 1: (a), (b), (c) The KEK data on the reaction tt+tt" tt^tt" [3]. (a) The normahzed S wave 
intensity l^sp. (b) The normahzed D wave intensity |>i_Dp; the curve is the fit using Eq. (1) with the 
parameters of the /2(1270) presented in Eq. (2). (c) The relative phase S between the amphtudes As 
and Ao. (d) The I = mr S wave phase shift Sq obtained from the data on \As\'^ alone under the 
assumption tjq is unity. 
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Figure 2: The I = 2 irir S wave phase shift Sq. The data are from Refs. [17] (soUd circles) and [18] 
(open circles). The curve corresponds to the fit described in the text. 
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Figure 3: The solid circles show the values of the phase sift Sq (a) and inelasticity tjq (b) extracted 
from the KEK data [3] on the reaction tt+tt" — s- tt^tt^ in the case that the m dependence of the phase S 
corresponds to variant (i) described in the text, (c), (d) The same for variant (ii). The open diamonds 
show the values of Sq corresponding to Fig. Id. 
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Figure 4: The results of the extrapolation of the BNL data [4]. (a) The extrapolated and normalized 
S wave intensity, (c) The extrapolated and normalized D wave intensity, (e) The extrapolated relative 
phase 6 between the S and D wave amplitudes. The slopes bs (b), bu^ (d), and a (f) as functions of 
m. The solid circles correspond to the physical solution. The open circles [and also the open triangles in 
plot (e) for S] correspond to the other solution. The lower and upper curves in plot (c) are the fits using 
Eq. (1) with the parameters of the /2(1270) presented in Eqs. (8) and (9), respectively. 
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Figure 5: The phase shift 6q and inelasticity r/g extracted from the BNL data [4]. Plots (a) and (b) 

correspond to the physical solution (see Fig. 4) with 5 > for all to. Plots (c) and (d) correspond to the 
physical solution with a transition of the phase (5 at m « 1 GeV from the branch pertaining to its positive 
values to that with its negative ones (see Fig. 4e). Plots (e) and (f) correspond to the combination of 
the physical solution with (5 > for to < 1 GcV and the other solution for to > 1 GeV with (5 > in 
the region 1 < to < 1.28 GeV and with (5 < in the region 1.28 < to < 1.52 GcV (see Figs. 4a and 4e). 
Plots (g) and (h) correspond to the combination of the physical solution with 5 > for to < 1 GeV and 
the other solution for to > 1 GeV with (5 < in the region 1 < to < 1.52 GcV (sec Figs. 4a and 4e). The 
open circles show the values of the phase shift Sq obtained from the data on \As\^ for to < 1 GeV (see 
Fig. 4a) alone under the assumption ryg is unity. 
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